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ABSTRACT 

We determine the NLO chiral low-energy constant L\q, and combinations 
CI 2 C'ei + ^805 ^13 ~ ^62 + ^ 8 i! ^ 61 ; ^^^1 Cgy, of the NNLO chiral low-energy 
constants incorporating recently revised ALEPH results for the non-strange vec¬ 
tor [V) and axial-vector (A) hadronic r decay distributions and recently updated 
RBC/UKQCD lattice data for the non-strange V — A two-point function. In the 
MS scheme, at renormalization scale fi = 770 MeV, we hnd L\q = —0.00350(17), 
C[2 + Ch + Qo = 0.00237(16) GeV-^ - C^i + Qo = -0.00056(15) GeV^ 
Cn - CI 2 + Qi = 0.00046(9) GeV^ = 0.00146(15) GeV^ and Cl^ = 
0.00510(22) GeV~^. With errors here at or below the level expected for con¬ 
tributions of yet higher order in the chiral expansion, the analysis exhausts the 
possibilities of what can be meaningfully achieved in an NNLO analysis. We 
also consider the dimension six and dimension eight coefficients in the operator 
product expansion in the V — A channel. 


I. INTRODUCTION 


Recently, a revised version of the ALEPH data for the non-strange vector {V) and 
axial vector (A) spectral distribntions obtained from measurements of hadronic r decays 
became available. These results corrected a problem, uncovered in Ref. j^, in the publicR 
posted 2005 and 2008 versions of the correlations between different energy bins.^ In Ref. 

’), as 


isim: 

li- 


This 


we analyzed these data in order to extract the strong coupling at the r mass, 
well as OPE condensates, following the strategy previously developed in Refs, 
analysis leads to a complete description of the V and A spectral functions as a function of 
the energy-squared s, including the region s > rri^. 

Complete knowledge of the spectral functions allows one to construct a representation 
of the (subtracted) vacuum polarizations in both channels, and the unsubtracted vacuum 
polarization in the V — A channel, as a function of s. Combining these representations with 
the analytic expressions derived from chira]_j)erturbation theory (ChPT), which are known 


a^6 

to next-to-next-to-leading order (NNLO) [7| allows us to extract the low-energy constant 
(EEC) Cgy, and a linear combination of Lio, C 12 — Cqi -|- Cgo and C 13 — Cq 2 + C'si from 
conventional chiral sum rules for the non-strange V — A channel @, [^, and C 12 + C^i + Cgo 
and a linear combination of Lio and C 12 — Cqi + Cgo from flavor-breaking {ud — us) chiral 


sum rules 
of L 5 and Lg 


10| in the V ± A channels. These determinations employ as input existing values 


)loy 

I3i 


estimates existing from both phenomenology I 111 . 1121 
In order to disentangle further Lio, C 12 — Cqi -|- Cgo 


and the lattice [13|, [14 
and Ci 3 — C 62 + Cgi, we exploit the 
dependence of the coefficients of these EEC combinations appearing in the V—A polarization 


on the pion and kaon masses using data for this polarization from the lattice |l5l. II6 

The goal of this article is to update the analysis of Refs. [ii®, replacing the experi¬ 


mental data for the non-strange s 
with the revised data from Ref. 
Ref. 


16 


spec 

• a, 


ectral functions, which previously came from OPAE [18 


m. 


15 


from 
In 


and, at the same time, updating the lattice input of 
The expectation is that the errors on Lio and the accessible NNEO combinations 
will decrease, because the AEEPH data are more precise, especially in the low-energy region. 
Improvements in the lattice data should also help in the process of disentangling L 
the combinations C 12 — Cei -|- Cgo and C 13 — C 62 + Cgi, described previously in Ref. 
the current analysis, we emplw a slightly different input for L 5 , choosing now to use the 
2 -I- l-ffavor estimate of Ref. [l^, LUu = 770 MeV) = 0.84(38) x 10“^. This value is the 
2 -|- 1-flavor estimate adopted in Ref. 1^ , and straddles several nominally more precise, but 
not mutually consistent, estimates, including the result = 770 MeV) = 0.58(13) x 10“^ 
of Ref. im used in our previous analysis. 

While the main emphasis here is on the EECs of ChPT, we will also update our estimates 
for the operator product expansion (OPE) condensates C^y-A and Cgy-^, which are order 
parameters (in contrast to the analogous condensates in the V + A channel). By their nature, 
these condensates are sensitive to the high-energy part of the spectral function, which is less 
well known. Therefore, it is interesting to compare the values for the condensates obtained 
from the OPAE and AEEPH data. 
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The corrected data may be found at http://aleph.web.lal.in2p3.fr/tau/specfunl3.html. 

We hope to revisit the determination of L 5 and Lg from the ud and us V polarizations in the future, but 
it remains to be seen whether the errors from such an analysis are competitive with those of Refs. Illll4|. 


Here we will use values from the literature. 
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II. THEORY COMPENDIUM 


In this section, we briefly collect the deflnitions and relations between quantities needed 
in order to present our results. More detailed overviews can be found in Refs. 1^, [iTj] and 
references therein. We first consider the required sum-rule results, and then collect relevant 
results from ChPT. 


A. Sum rules 


The vacuum polarizations, (Q^), with T = V, V it A and J = 0,1, are the spin 
J scalar parts of the corresponding current-current two-point functions, and are 

related to the corresponding spectral functions, p^\s), by finite-energy sum rules (FESRs). 
The flavor ud and us versions of these spectral functions can be obtained from experimental 
differential hadronic r decay distribution data for energies up to the r mass. Above the r 
mass one needs a theoretical representation, and we will use the one obtained in Ref. j^. 

First, let us consider the non-strange V — A channel, for which the vacuum polarization 
obeys the unsubtracted dispersion relation 

ny_A(Q^) = nK>(0t = j"ds ~ ■ (2.1) 

~ (9'“' - . 

where the Euclidean momentum-squared = —q^. Here pv (pa) is the non-strange, 
J = 1 vector (axial) spectral function summing the angular momentum J = 1 and J = 0 
contributions. Generalizing the definition of we also define functions fo be 

used in the restricted sense employed below, involving additional polynomial weight factors 
w{y)\ 

nlrh(Qt = His w(sla„) (2.2) 

Jo s + V 

with 0 < So < In what follows, we will use the weights 

Wk{y) = {l-yf, /c = l,2. (2.3) 


In evaluating the integral in Eq. fl2.2p . we will use the ALEPH experimental data for the 
spectral functions for s < sq, and the duality-violating (DV) part 

Pv{s) - pa{s) ^ Py^(s) - p5^(s) , s > So , (2.4) 


above so, with Py^(s) — Pa^(s) equal to I/tt times the imaginary part of defined, 

in turn, from 

t 1 v-a(Q^) = n°™ (Q'l) + . (2.6) 

for \Q‘^\ > So, where the OPE part has the form 


k=l 


C2k,VIA 


( 2 . 6 ) 
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with C 2 k,viA the OPE coefficients. We will always assume that we can choose sq smaller 
than but large enough that the separation fl2.5l) into OPE and DV parts makes sense. 
We use a model for the DV parts of the spectral functions: 


Pv/a{s) = e ^ sin (a +/9s), 


(2.7) 


with a, /9, 7 , and 6 parameters which differ in the V and A channels. The form of this 
ansatz was motivated in Ref. 1^, and it was shown in Refs. j3, [sj that this model can be 
used to successfully parametrize the resonance structure in the data for si^l.4 GeV^. Here 
we will fix the values of the V and A DV parameters using the results of the FOPT fit of 
Table IV of Ref. with Smin = 1.55 GeV^.^ 

In what follows we will denote by H^, H^^ and the versions of Hyi, H^^ and pa 
from which the pion pole contribution has been subtracted. Analogous pion-pole-subtracted 
versions of the V i: A polarizations and spectral functions are denoted Ily±A = Ify i If^, 

^v±A = nV ' ± and = pv± Pa- 

For < 4m^, IIv-a{Q^) admits a Taylor expansion, and we can thus dehne the intercept 
and slope at = 0, 


Uv-a{Q‘') = -SLfo - ISGstQ" + 0{Q^) 


( 2 . 8 ) 


Employing FESRs for the weights fITSp and analytic results for the OPE coefficients C 2 y-A 
and Cyv-A 21|, it follows that 


-sLf, = nv_^(o) = nir^J,(o) + 


VI 

So 


(2,9) 


= n7>) + k£ 

■So 


17 f as{so) V mu{so)md{so) _ ^ A 4 q;^(so) 

levr^ V tt / /2 2so V 3 


TT 


Since the terms proportional to q;s(so) lead to effects smaller than the experimental errors, 
we will omit these terms from the actual analysis leading to our results for and . 

In addition to the information obtained from the flavor ud V-A channel, further con¬ 
straints can be obtained from inverse-moment FESRs (IMFESRs) for the flavor-breaking 
differences 


AnT(Q^) = ni“7>(Q" 


n, 


( 0 + 1 ) 

us\T 


(Q^ 


( 2 , 10 ) 


defined in Ref. [17|. Note that the An/r(Q^) are finite and satisfy unsubtracted dispersion 


relations. The IMFESRs of Ref. [17] have the forms 

w{s/so) . , . 1 


'■so 


Anv(o) = 

Any±yi(o) = 


ds ■ 


■ Apy(s) -h 


' sth 

P w(s/so) , , , 

/ ds -Apy±^(s)± 

Gth 

1 f . w{z/so) 


271 i 


|2|=S0 

k 


dz'AlM. 


— frJvK) - -/”(!/.) 

So So 


J TT 


2771 


dz 


|2|=S0 


[Any±^(—2;)] 


OPE 


( 2 . 11 ) 


Essentially identical results are obtained if one employs instead the results of the corresponding GIFT fit. 
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in which Sth is the continuum threshold 4m^, = ml/so, Hk = ''^Ir/so, and 

f{y) = ‘^{w{^)-w{y)) . (2.12) 

As long as we retain the exact An(— z) in the contour integrals of Eq. (12.lip , the full right- 
hand sides are necessarily independent of the weight choice tc, provided we restrict our 
attention to w{y) all having t(;(0) = 1. In Eq. fl2.1ip we dropped the DV term from the 
split fl2.5p . and kept only the OPE contribution. It is reasonable to do so, because the only 
weights we will use will be triply pinched, he., contain a factor {z — Sq)^, which suppresses 
DVs strongly. In our analysis, we will use the weights 


w{y) = {1-yf , (2.13) 

WDK{y) = (1 - 2/)^ (^1 + 2/ + ^2/^^ =l-2y+^{y^+ y^+ y^-y^) . 


The IMFESR with tcoK was first considered in Ref. 10 


The quantities ny_yi(0), ny_^(0), Anv/(0) and Any±^(0) can be obtained from hadronic 
r-decay data, which yield the spectral functions pv/a{s), either directly, for s < m^, or, where 
it is needed for s > in the V — A channel analysis, indirectly through Eq. fl2.4p . using 


the values of the DV parameters obtained from the hts to these same data, described in 
Ref. 01 . In addition, one needs the experimental values of tUtt, fn,^ mR and fK- For a 
detailed discussion of the OPE contributions to Eq. fl2.1ip . we refer to Ref. |l7|.® A key 
point is that the numerical contributions from the OPE terms to Anv/(0) and Any±yi(0) 
are very small. That implies that even if these OPE contributions are not very well known, 
and one has therefore to include very conservative estimates of their errors in the total error 
budget, the impact on our final errors is minor. 

As already noted above, we have dropped DV contributions in Eq. fl2.11l) . The reason is 
that these are very suppressed for the weights fl2.13l) . which are triply pinched, and moreover 
suppress the large-s region by an additional factor 1/s. Since the Sq dependence from both 
the OPE and DVs is non-trivial, our treatment of the OPE and the omission of DVs can 
be tested by checking the Sq independence of Any(O), AnT/+yi(0) and Any_A(0) produced 
using the right-hand side of the corresponding IMFESR. 


B. Chiral perturbation theory 

The motivation for considering the quantities L®q, OI®, Any(O), and Any±A(0) is that 
they all depend on NLO and NNLO LECs of the chiral effective theory, and thus yield 
information on the QCD values for these LECs if sufficiently accurate data (from experiment 
or lattice QCD) are available. Here we will collect the relevant NNLO ChPT expressions 
needed in order to connect the quantities defined in the previous subsection to the LECs of 
ChPT. 


^ Our normalization is such that /jr = 92.21(14) MeV. 
® See, in particular, Sec. III.B of that reference. 
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(2.14) 


The representation of IIv-a{Q^) to NNLO in ChPT has the form 


Hv-Am = + + + 

+32(2m^ + ml) (^[3 - + C;,) - IGC^^tQ" + ^ 9 ) , 

where the explicit expression for Tg) can be reconstitnted from the resnlts of 

Ref. The snbscript vriP indicates that R depends also on mK and /^r, in addition 
to the explicitly shown argnments.^ L®® ^^^1 are then given by: 


nv-^(O) = —8L 


eff 

10 


~ 4(2/i7r + + 16(2/i^ + ^K)Lg 

+32ml {Cl, - Q, + C;,) + 32{2ml + ml) {Cl, - C;, + C;,) 

-\-R.j^k{^^i 0 ) 5 


ny_A(0) = 


■lec'st 


- 16 Q 7 + 
, dRT,K 


47r2/2 


Li^ 1 

1 - log — + - log 
mi 3 


m 


K 


mt 


ml 


Pp — 


327r2/i 




log 


Q2=0 


mt 




(2.15a) 


(2.15b) 


(2.15c) 


where RT,K{f^,Q‘^) is the part of RTrK{f^,Q‘^', L^) independent of Lg. Here the snperscript r 
denotes the valnes of LECs renormalized in the MS scheme at scale p, which we will take 
to be /i = 770 MeV in what follows. 

We will also need the NNLO ChPT expressions for Anp(O), T = R, R ± H, bnt only 
at physical valnes of m 7 r, m^ and f.^. We therefore give the expressions in terms of the 
LECs with the nnmerical valnes of the coefficients for the NLO LECs g evalnated at 
m^ = 139.57 MeV, Mk = 495.65 MeV,® U = 92.21 MeV, and p = 770 MeV: 

Any(O) = 0.00775 - 0.7218 + 1.423 + 1.062 L^g + 32(m^ - m^)^^^ , 

Any+A(0) = 0.00880 - 0.7218 + 1A23 + 32{m\ - ml) [Cl, + C^^ + C'g'o] , 

Any_A(0) = 0.00670 - 0.7218 + 1.423 L(; + 2.125 L^g 

~32{mj^ — m^)[Ci, —Cqi +C,(^ . (2-16) 


Eqnations fl2.14p . fl2.15p and fl2.16p give access to combinations involving the LECs L^g, 
Cgy, and the linear combinations 

^ 32m2 (C[2 - + Qg) , 

Cl = 32 (ml + 2ml.) (CI 3 - Cl, + CJ,) , (2.17) 


® Since the value of Lg is well known 


12| , we treat the loop contribution proportional to this LEG as known, 
and we thus include this contribution in (m, Lg). 

In Ref. 


15| we denoted this term simply as R{Q^i Lg). 


This is the average of the charged and neutral kaon masses. Taking just the charged or neutral value has 


no impact on our final results. 
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In addition, the LECs Lg, Cg^ 


and the linear combination C'i 2 +C' 6 i +^80 appear in Eq. fl2.16p . 
Onr goal is to use the ALEPH hadronic r-decay data to extract L^q, Cg^^, C^j and Cl 2 ^Cli + 
(PgQ, using the known values of Lg and Lg. However, L^q appears in a linear combination 
with Cq and C[ that does not depend on and we thus need other information in order 
to disentangle L^q, Cq and C[ from each other. This can be done by using lattice data for 
different values of the meson masses, and such data are available for T = V — A [i3.® 

We thus will consider, following Ref. 15 




(2,18) 


the difference between the non-strange pion-pole-subtracted V — A correlator 
evaluated on a lattice ensemble B with tt and K masses and decay constants different from 
the physical ones, and the same correlator, evaluated at the same Q^, for the 

physical quark mass case. The latter is obtained from the dispersive representation using 
spectral functions obtained from the hadronic r decay data. It then follows that in terms of 
LECs this difference can be expressed as 




L,E 


rr 

-^10 


+ 


rL,E^r 
^0 ‘-'0 


+ s';'^ci 


(2,19) 


where Lg) and the Q^-independent coefficients are known in terms of the 

lattice and physical meson masses and decay constants, and the chiral renormalization scale 
/i. Evaluating Eq. fl2.18p using lattice and dispersive results, Eq. fl2.19p yields a constraint 
on L^q, Cq and C[ for each ensemble, E, and each Q^. As explained in more detail below, 
different choices of for fixed E provide self-consistency checks on the use of the lattice 
data. 


III. INPUT DATA 


We will evaluate and using ALEPH experimental data [H for the 

spectral functions py{s) and p^(s) for s < sq = Sswitch/^ and approximating the difference 
py{s) — Pa{s) by Eq. fl2.4p for s > Sq =^witch 5 with values for the DV parameters from the 
combined V and A channel fits of Ref. . We will choose Sswitch to be the upper end of an 
ALEPH bin, obtaining, in the notation adopted in Ref. Q, 


f(«i) 


KUQ') = 


/ sbin+dsbin/2 

ds w 

bin-dsbin/2 


py(sbin) — p^(sbin) 


+ 


ds w 


•^switch 


■^switch 


Sswitch/ sbin + Q2 

Pr(s)-pS^(s) 

s + g2 


(3.1) 


Here we will choose Sswitch = 1-55 GeV^, the value of Smin which produced the best £t to the 
weighted spectral integrals in our extraction of as in Ref. [3]. Since we are only interested in 

the behavior of A^ylj^{Q^) at values of -C 1.55 GeV^, the right-hand side of Eq. fld.ip is 


® No such lattice analysis has been performed yet for the us channel, which is why we will use known values 
for Lg and Lg here. 

For details on the handling of the ALEPH data, see Ref. [^ . 
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very insensitive to the precise choice of ^witch, and varying it within the range of Smin values 
for which we obtained good fits in Ref. has no effect on either the values or errors we will 
obtain for the LECs below. 

We have fully propagated all errors and correlations in the results we will report on below. 
In particular, the DV parameter values used in Eq. fl3.ip are correlated with the data, and 
we have computed these correlations using the linear error propagation method summarized 
in the appendix of Ref. (see, in particular, Eq. (A.4) of that reference, which can be used 
to express the parameter-data covariances in terms of the data covariance matrix). 

For the us data needed in order to evaluate the strange spectral integrals in Eq. fl2.1ip we 
will use exactly the same treatment and input as used in Ref. [l^, and we refer to Sec. III.C 
of that article for details. We collect the values of all other external input parameters: 


rriK 

rrin 

U 

Ik 

Ll{fi = 770 MeV) 
Ll{jj = 770 MeV) 


139.57 MeV , 
495.65 MeV , 
547.85 MeV , 
92.21(14) MeV , 
110.5(5) MeV , 
0.84(38) X 10-3 , 
5.93(43) X 10-3 , 



(3.2) 


For the value of we choose the value reported of Ref. [l^, which has a larger error than the 
value of Ref. H we used in Ref. [l^. The comparison with other values in Ref. [l3[ shows 
that the error quoted in Eq. 03.21) above covers these other values, and we thus consider its 
use to represent a conservative choice. 


IV. RESULTS 

In this section, we present the results of our analysis, dividing the presentation into several 
parts. We first present our values for L®g and C*!®, which are based purely on non-strange 
r-decay data, then derive additional constraints employing also the lattice data, and finally 
use the us spectral data to obtain further constraints via the flavor-breaking IMFESRs. The 
output is a determination of L[q, Cq^, Cgj, CI 2 — C^i + C'so ^i3 ~ ^62 + C'sii in terms of 
the known values for Lg and Lg. In the final subsection, we switch gears, and consider what 
we can learn from the ALEPH data about the OPE coefficients C^y^A and C^y^A- 


A. Effective LECs 


We hrst give the values of L®q and Ol®, which follow directly from Eq. fl2.9p and the 
evaluation of Eq. fl3.1l) . For the three different weights we hnd 


' efF 
^10 


6.482(64) X 10-3 , 

(from ny_A(0)) , 

(4.1a) 

6.486(64) X 10-3 , 

(from ny_)Ji(0)) , 

(4.1b) 

6.446(50) X 10-3 , 

(from ni;'_!li(0)) . 

(4.1c) 


These values are consistent with those found employing OPAL data in Ref. [15[, but more 


precise, with errors about a factor two smaller. The contribution from the DV part of the 























spectral function in Eq. fl3.ip ranges from 3% for the first estimate in Eq. fl4.ip to about 1% 
for the third estimate. Their size is thus comparable with the quoted errors, suggesting that 
the uncertainty in the DV part due to the use of a model for DVs, Eq. fl2.7p . is negligible. 
From the slope of Yiv-AiQ"^) at zero we find 

= 8.38(18) X 10“^ GeV^ . (4.2) 


In this case, the contribution from the DV part is about 1%. These results are in good 
agreement with those in Refs. 23|. However, our errors are somewhat smaller, and 
our analysis employs versions of the DV contributions fitted individually in the V and A 
channels, as required by the data, avoiding the simplifications employed in Ref. 22|, l23(] (see 


also the discussion in Ref. |15|). 


From our best value for L®q, Eq. fl4.1cl) . and using Eq. fl2.15al) . we find the constraint 

= 0.6576L^{^ + 0.001161 - 0.1712L(; + 0.08220(C(; + C[) (4.3) 

= -0.004094(33)ALEPH(74)i. + 0.08220(Co^ + C[) . 

Together with information from the lattice on the combinations Cg and C[, this constraint 
will yield an estimate for L^g. Likewise, Eq. (I2.15bp leads to the constraint [HI 


^eff _ 

^87 “ ^87 


and, with Eqs. (14. 2 p and 


0.292L^ + 0.00155 GeV 
, this leads to the estimate 
= 5.10(22) X 10-3 GeV-2 


-2 


(4.4) 


(4.5) 


B. Constraints using the lattice 


As shown in Ref. [l^, and noted already above, useful independent constraints on L^g, 
Cq and C[ can be obtained by considering the difference of the lattice-ensemble, E, and 
physical-mass, continuum results for IIv-a{Q^), evaluated at the same Q^. In what follows, 
we will use the superscript phys to specify quantities evaluated with physical values of the 
masses and decay constants. 

Recasting the NNLO representation, Eq. fl2.19p . in the form 


rL,E rr I 


cL,E/-tr 
^0 ^0 


+ = 


Ant’!^(g2) 


(4.6) 


one notes that the left-hand side is explicitly Q^-independent, while the right-hand side is the 
difference of two Q^-dependent terms. For values of for which the NNLO representation 
is reliable, the for hxed ensemble E but different should be compatible within 

errors, thus providing non-trivial self-consistency checks. We will denote the average of the 
for an ensemble satisfying these self-consistency tests by 
The explicit expression for AR^’^{p, Lg) = ~ Q^; Lg) follows 

from the results of Ref. [^, but is very lengthy and hence not given here. The explicit 
expressions for the mass-dependent constants appearing on the LHS of Eq. fl4.6p are 


5 

(5 

(5 


L,E 

10 ~ 

32 (2p, + - 

L,E _ 

K]LE 

0 — 

[ml]P^y^ ’ 

L,E _ 

[ml + 2mlt’^ ^ 

1 — 

[ml + 2m‘j^f’^^^ 


(2/i7r “|- Pk') 


phys 




(4.7) 
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The lattice data we employ in forming the lattice-minus-continuum constraints are ob¬ 
tained using the nj = 2 + 1 domain wall fermion ensembles of the RBC/UKQCD collabo¬ 
ration. Details of the underlying simulations may be found in Refs. [ 2 ^ 25|, with updated 
values of the lattice spacings a, obtained after incorporating results from the new physical 
point ensembles, given in Ref. j26| . 

We have used the following criteria in deciding on the choice of ensembles and values to 
be employed. First, we restrict our attention to ensembles with m.^ < 350 MeV. Second, we 
require the ensemble to have sufficiently many points in the expected range of validity of 
the representation, Eq. fl2.19p . that meaningful self-consistency tests can be performed. With 
the errors at the lowest point turning out to be very large for all ensembles considered, 
this means that a minimum of two additional such Q^, or three in total, are required. 
Finally, we identify the range of validity of the representation Eq. fl2.19p as follows. We first 
note that the supplemented NNLO representation of liy -discussed in more detail in 
Refs. HQ, and the corresponding NNLO representation given above, both yield the same 
representation of the lattice-continuum difference The supplemented form of 

IIv-a{Q^) incorporates resonance-induced NNNLO contributions analogous to those already 
encoded in the NNLO contribution proportional to through the inclusion of an additional 
analytic term The inclusion of such a term was shown to extend the reliability of the 

supplemented version of the representations to significantly larger in both the V — A 


and V correlator cases [15|, l27|. Here we investigate the supplemented NNLO fit by fixing 


and Cgf to the values given in Eqs. fl4.ip and fl4.2p . and htting the additional effective 
mass-independent NNNLO LEG, C, to the dispersive results for ny-yi(Q^) in the window 
0 < < Qmax- The range of validity of the supplemented form is then identihed as the 

largest which such a £t is successful within errors. The results of this exploration 

show that the supplemented form can be employed up to 0.25 GeV^, but not beyond. We 
thus restrict our attention further to ensembles having at least three in the region below 
0.25 GeVl 

Three RBG/UKQGD ensembles satisfy the above criteria, two coarse ensembles with 
1/a = 1.379(7) GeV and 171 -!^ = 172 and 250 MeV, and one intermediate ensemble with 
1/a = 1.785(5) GeV and = 340 MeV jl^. These are labelled E = 1, 2 and 3 in what 
follows. The coarse ensembles have eight points each below 0.25 GeV^, the intermediate 
ensemble has four. The results for the V — A correlators for all three of these ensembles 
pass the self-consistency tests discussed above. 

The lattice-continuum constraints for the first two cases were obtained previously in 
Ref. 1^. While the results of Ref. 2g lead to a small shift in the value of 1/a for these 
ensembles, this change affects the constraints only through the values of the at which 
the dispersive representation of the physical-mass correlator is evaluated, the values of 
being fixed by the lattice size. The small resulting changes turn out to have no impact on 
the resulting averages for these ensembles to the number of signihcant figures previously 
quoted. The resulting average are thus those given in Ref. 16l |. 


= 0.0007(17) , 
TL’2 ^ 0.0039(21) . 


(4.8) 


The third ensemble has a relative error on af-jr a factor of 2 smaller than that for the two 
coarse ensembles, and hence a smaller uncertainty in the pion-pole subtraction involved in 
forming from the directly measured unsubtracted version. In order to improve 


10 






















further the associated constraint, the statistics for this ensemble were increased using mul¬ 
tiple time sources. This increase in statistics was greatly aided by the use of the HDCG 
algorithm of Ref. , employed in performing the propagator inversions. The covariances of 
the corresponding lattice-minus-continunm differences for different (eqnal to the snm of 
the covariances of the corresponding lattice and dispersive results) are strongly dominated 
by the lattice contribntions. With the covariance matrix available, the average constraint 
value, T^’^, can be obtained by a standard, fully correlated ht to the four with 

< 0.25 GeV^ available for this ensemble. The resnlt is 

= 0.00625(48) . (4.9) 

The error here reflects only the errors on (and correlations among) the lattice results at 
different and dispersive results at different Q^. Additional errors dne to the uncertainties 
on the input quantities L^fi) and which are 100% correlated with the analogous 

uncertainties entering the ny_^(0) and flavor-breaking IMFESR constraints, are handled 
separately below. 

C. Constraints using IMFESRs 

In this subsection, we evaluate Anr(O) with T = V, and V ± A from Eq. fl2.11l) . for 
values of sq between 2.15 GeV^ and m^. Of course, AnT(O) is independent of Sq, implying 
that evalnating the right-hand side of the expressions in Eq. fl2.1ip and checking for sq 
independence provides a self-consistency check on the use of the OPE, and the assnmption 
that DVs can be neglected. The valnes we hnd are, using w(y) = WDK(y) in Eq. (12.lip . 

Any(O) = 0.0224(9) , (4.10) 

AnA(O) = 0.0113(8) , 

Any+A(0) = 0.0338(10) , 

Any_A(o) = 0.0111(11). 

The result for An^(O) obtained from the analogous A channel IMFESR is also included 
for completeness. While the axial case is not independent of the others, performing the 
axial analysis directly is the most straightforward way to take into account the correlations 
amongst the other channels and arrive at the correct error for the A case. The quoted 
errors take into account the experimental errors in the ALEPH data, the uncertainties in 
the estimates of the OPE contribution, and the small residual So-dependence observed as Sq 
is varied over the analysis window noted above. As a further check of the self-consistency of 
the valnes in Eq. (I4.10p . we have rerun the analysis using w(y) = w(y) in place of WDK(y), 
and hnd results compatible with those obtained using WDK(y) to well within the errors 
qnoted in Eq. fl4.10p . The analysis method leading to the valnes fl4.10p is identical to that 
of Ref. [Isl , to which we refer for a detailed discnssion. 

From the valne for Any+^(0) and Eq. fl2.16p we hnd, using the valnes of and Lg in 
Eq. 03.2p . 

+ C^so = 0.00237(13)^n(4)Ls(8)L5 GeV^ (4.11) 

= 0.00237(16) GeV^ , 

where the subscript AH refers to the error coming from the value for Any_|_^(0) in Eq. 04.101) . 
and on the second line we have combined errors in qnadrature. The sum rules for Any_^(0) 
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and Anv^(O) provide two independent constraints involving combinations of L\q and other 
NNLO LECs: 


2.125I/;^q 32(m^ m^){Ci2 C'g;^ + Cgg) 


and 


-0.0034(ll)^n(3)Le(6)L5 GeV-2 
-0.0034(13) GeV^ , (4.12) 


1.062L^o + 32(m^ - m^)^, = 0.0071(9)^n(3)Ls(6)L5 GeV" (4.13) 

= 0.0068(11) GeV^ . 

Employing the constraints Eqs. fl4.3p . fl4.12p and the E = 1,2,3 versions of Eq. fl4.6p . 
with Eqs. fl4.8p and fl4.9p as inpnt on the right-hand side and Eq. fl4.7p as inpnt on the 
left-hand side, we hnd 

L\^ = -0.00350(ll)(13)L5,Lg = -0.00350(17) , (4.14) 

Cl = -0.00035(9)(4)Lg,Lg = -0.00035(10) , 

Cl = 0.0075(13)(8)l5,L9 = 0.0075(15) . 


Here the first error quoted comes primarily from the errors on the lattice data,^^ while the 
second error comes from the errors in Lg and Lg. The value of Cq translates directly into 



CI 2 - Ggi + Clo = -0.00056(15) GeV^ , 

(4.15) 

the value of C[ into 

CI 3 - CI 2 + = 0.00046(9) GeV^ , 

(4.16) 

while the value of L[q 

in Eq. f|4.14p toeether with Eq. f|4.13p leads to the estimate 



di = 0.00146(15) GeV^ . 

(4.17) 

For completeness, the result for the NNLO LEG combination entering the flavor-breaking A 
IMFESR is 

+ do = 0.00090(9) GeV^ . (4.18) 


D. V — A condensates 

In addition to the LECs extracted in the previous subsections, Y\.v-a{Q‘^) also provides 
information on the OPE coefficients C^y^A and Cgy-A- Ad^ting Eq. (4.20) of Ref. [isj to 
the case of the ALEPH data, following the notation of Ref. used also in Eq. fl3.ip above. 


Lattice errors in the differences (14.101) dominate over the errors coming from the r spectral data. 
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these two coefficients are given by 


C(iy-A — 

sbin<Ssw 

-2fl [ml - Sswitch)^ + 
Cs,V-A = — ^ f 


/•sbinH-dsbiii/2 

' ds {s- Sswitch)^ (pi/(sbin) - p^(sbin)) 

sbin—dsbin/2 


(4.19a) 


' ds {s- Sswitch)^ (Py^('S) - Pa^('S)) , 

^switch 

sbinH-dsbiii/2 

_ ds {s- Ss„itch)^(s + 2Sswitch) (py(sbin) - p^(sbin)) 

sbin<s,„it,h 

(^TT "Sswitch) 4” 


ds {s - Sswitch)^(s + 2Ss„itch) {pv^{s) - p5^(s)) 


(4.19b) 


•^switch 


The first of these two expressions involves contributions proportional to C 2 k,v-A for k = 
1, 2, 3, but the leading-order expressions [1^, 2l[ 


C2,V-A- 


f 2\ f 2 

-— m„(/i )md[p 


TT"* 


) 1 




TT 


17 

T 


log — + c 

II 


+ 


_ 8 2 2 I 

Ciy-A——-;: —+ 
3 71 


(4.20a) 

(4.20b) 


with c a numerical constant whose value is not required in what follows, suggest that the 
D = 2 and D = A terms are numerically negligible. A similar observation holds for the 
second of these expressions. The use of Eq. fl4.19l) to estimate C^y^A and C^y-A thus 
rests on the assumption that these estimates for the D = 2 and D = A contributions are 
sufficiently reliable, which is equivalent to assuming that the two Weinberg sum rules hold 
exactly. In this case, Eq. fl4.19p leads to the estimates 


ALEPH: 


Cey_A = (-3.16 ±0.91) X 10-^ GeV® 
Csy-A = (-13.0 ± 5.5) X 10"^ GeV® . 


(4.21) 


These results correspond to the choice Sgwitch = 2.2 GeV^, which yields the smallest estimate 
for the errors on C^y^A and Cgy^A- We have checked that the central values remain stable 
as a function of Sswitch as Ss„itch is varied between Smin = 1.55 GeV^ and ml. The results of 


Eq. 04.21 p are in agreement within errors with those of Refs. j22|, [23 


Let us compare these values with those we found from the OPAL data in Ref. [15 

OPAL : Cey-A = (-6.6 ± 1.1) x 10“^ GeV® , 

Csy-A = (5 ± 5) X 10-3 GeV® . 


(4.22) 


These values differ by 2.4 a from the ALEPH values in Eq. 04.2ip . 

Instead of the weights employed in Eq. 04.19p . one can use the weight to estimate 
C^y-A, and the weight to estimate Csy-A- If we do so, we find the values 


Cey-A = (-4.4 ± 5.2) x lO’^ GeV® 
Csy-A = (-8 ± 18) X 10-3 GeV® . 


(4.23) 


This suggests that central values of C 2 y-A and C 4 ^y-A produced by the ALEPH data are 
larger than those following from Eq. 04.201) . Direct determination of these coefficients from 
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the data, using the weights 1 and s on the right-hand side of Eq. fl4.19p . yield values C 2 y-A = 
(0.14 ± 0.24) X 10“^ GeV^ and C^y-A = (0.1 ± 1.2) x 10“^ GeV^. These central values are 
indeed larger and explain the difference between the estimates fl4.19p and fl4.23p . though 
the results are consistent with Eq. fl4.20p within errors. The smaller errors reported in 
Eq. fl4.22p are thus a consequence of the assumption that C 2 y-A and G 4 y can be neglected, 
or, equivalently, the assumption that the first and second Weinberg sum rules are exactly 
satished. 

In view of the discussion above, we conclude that current data lead to somewhat con¬ 
flicting estimates for C^y-A and Csy-A- This is not surprising, because these coefficients 
are sensitive to the large-s region of the data. In addition, we observe the contribution from 
the DV terms to the expressions on the right-hand side of Eq. fl4.19p . he., the size of the 
DV contributions which must be subtracted to arrive at the values reported in Eq. fl4.2ip . 
are non-negligible: about —0.28 x 10“^ GeV® and 3 x 10“^ GeV® for C^y^A and Ggy-A, 
respectively. 


V. CONCLUSION 

Using the recently updated and corrected ALEPH results for the non-strange V and A 
hadronic r decay distributions, existing results for the corresponding strange decay distri¬ 
butions, and updated lattice results for the non-strange V — A polarization at heavier than 
physical meson masses, we have produced improved determinations of the NLO chiral LEG 
and a number of NNLO LEG combinations. Those results are given in Eqs. fl4.1ip and 
fl4.14p to fl4.18p . We have also used the non-strange ALEPH data to extract the dimension 
six and eight condensates appearing in the OPE representation of the V — A polarization. 
These results are given in Eq. 04.211) . 

The improvements produced by using the new lattice results and the new ALEPH data 
in place of the old OPAL data reduce the ht component of the errors on L\q and the 1/A'"c- 
suppressed NNLO LEG combination Gj'g — Cq 2 + Gg^ by a factor of roughly 2.5 compared to 
our earlier analysis. The fit errors on the remaining NNLO LEGs are about 2/3 of those of 


the previous analysis. Taking, as in Ref. [15[, 25% as an estimate of the expected reduction 


in size of contributions in going from one order to the next in the chiral expansion, we would 
expect the uncertainties from the neglect of NNNLO and higher contributions to be roughly 
6 % for and 25% for the NNLO LEGs. With the current fit errors, the total errors on 
all the LEGs determined have reached these levels, suggesting that the optimal practical 
precision one can expect for an NNLO analysis has now been attained. 

We hnd, in contrast, that using the higher precision ALEPH data produces essentially no 
improvement in the accuracy of the determination of the dimension six and eight V — A OPE 
condensates. Our results for these quantities are in agreement with those of other ALEPH- 
based analyses, and show about 2Aa discrepancies with the corresponding results obtained 
using the OPAL data. These discrepancies presumably reflect additional systematic uncer¬ 
tainties encountered in attempting to extract these small higher dimension contributions 
from existing data, and should be kept in mind if results based on one or the other of the 
two data sets are employed in other contexts. 
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